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' Abstract. We evaluate the matrix elements {OrP), where O — {1, /3, ianf]} are the standard Dnac 

I matrix operators and the angular brackets denote the quantum-mechanical average for the relativistic 

' Coulomb problem, in terms of generalized hypergeometric functions 3F2 (1) for all suitable powers. 

' Their connections with the Chebyshev and Hahn polynomials of a discrete variable are emphasized. 

As a result, we derive two sets of Pasternack-type matrix identities for these integrals, when p —^ 
—p — 1 and p — > —p — 3, respectively. Some applications to the theory of hydrogenlike relativistic 
• systems are reviewed. 

m ■ 

^! 1. Introduction 

a 

d ' Recent experimental progress has renewed interest in quantum electrodynamics of atomic hy- 
drogenlike systems. Experimentalists and theorists in atomic and particle physics are discovering 
I — \' problems of common interest with new ideas and methods. A current account of the status of this 



0^ 



fundamental area of quantum physics, which is more than a century old, is given in Refs. [29j, 
> ' [30], [SB], and [SH]- Exciting research topics vary from experimental testing of Quantum Electro- 
. dynamics (QED) to fruitful training models for the bound-state Quantum Chromodynamics and 
^ ', Bose-Einstein Condensation [H], [29], [30|, [3l], [38], [58], and [63] . 



The highly charged ions are an ideal testing ground for the strong-field bound-state QED. They 
• posses a strong static Coulomb field of the nucleus and a simple electronic structure which can be 
^ < accurately computed from first principles. It is possible nowadays to make massive highly charged 
ions with a strong nuclear charge and only one electron through the periodic table up to uranium. 



[ the most highly charged ion [26], [27]. These systems are truly relativistic and require the Dirac 
^ ' wave equation as a starting point in a detailed investigation of their spectra [38] , [55] . The binding 
. energy of a single /^-shell electron in the electric field of a uranium nucleus corresponds to roughly 

one third of the electron rest mass. For the simple hydrogen atom the nonrelativistic Schrodinger 

approximation can be used [S]. 

For the last decade, the two-time Green's function method of deriving formal expressions for the 
energy shift of a bound-state level of high-Z few-electron systems was developed [55] and numerical 
calculations of QED effects in heavy ions were performed with excellent agreement to current 
experimental data [26], [2^ (see [52], [53], [51], [57], [58], and [63] and references therein for more 
details). These advances motivate, among other technical things, evaluation of the expectation 
values (Or^) for the standard Dirac matrix operators O = {1, f3,ian(]} between the bound-state 
relativistic Coulomb wave functions. Special cases appear in calculations of the magnetic dipole 
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hyperfine splitting, the electric quadrupole hyperfine splitting, the anomalous Zeeman effect, and 
the relativistic recoil corrections in hydrogenlike ions (see, for example, [T], [5B], and [M] and 
references therein). We discuss convenient closed forms of these integrals in general and derive 
matrix symmetry relations among them which can be useful in the theory of relativistic Coulomb 
systems. 

The paper is organized as follows. In the next section we review the relativistic Coulomb wave 
functions and set up the notations. The expectation values (Or^) are evaluated in section 3 in 
terms of the generalized hypergeometric functions 3F2 (1) for all admissible powers of r. Their 
Pasternack-type matrix symmetry relations are established in section 4 and recurrence relations 
are given in section 5. We discuss special matrix elements and review some of their applications in 
the last section. An attempt to collect the available literature is made. The appendix A contains 
definition of the generalized hypergeometric series and proof of a required transformation identity. 
The appendix B deals with the Dirac matrices and inner product. 



2. Wave Functions for the Relativistic Coulomb Problem 



The exact solutions of the stationary Dirac equation 

Htfj = {cap + mc^l3 - Ze^/r)^ = Etfj (2.1) 

for the Coulomb potential can be obtained in the spherical coordinates. The energy levels were 
discovered in 1916 by Sommerfeld [59] from the "old" quantum theory and the corresponding 
(bispinor) Dirac wave functions were found later by Darwin [16] and Gordon [25] at the early age 
of discovery of the "new" wave mechanics (see also [9] for a modern discussion of "Sommerfeld's 
puzzle"). These classical results are nowadays included in all textbooks on relativistic quantum 
mechanics, quantum field theory and advanced texts on mathematical physics (see, for example, 
0, [3, [H], [SBI, [23, and [10]). The end result is 

F(r) 



X 



^yjm (n) 



G{r) 



(2.2) 



where the spinor spherical harmonics yf^ (n) = y. 



(i±i/2) 

' jm 



n) are given explicitly in terms of the 
ordinary spherical functions Yim (n) , n = n (6*, (^9) = r/r and the special Clebsch-Gordan coefficients 
with the spin 1/2 as follows [|, [7], [16], [6l]: 

/ 



yfm (n) 



'( j + l/2)T(m-l/2) 
2j + (1 ± 1) 

'(j + l/2)±(m + l/2) 



j±l/2, m-1/2 



n 



v 



2j 



:i±ii 



Y.. 



i±i/2, m+1/2 



n 



with the total angular momentum j = 1/2, 3/2, 5/2, ... and its projection m = —j, + 
(see also Section VI A of Ref. [60] for the properties of the spinor spherical harmonics). 



(2.3) 



The radial functions F (r) and G (r) can be presented as [lO 

Fir) \ _ a2/33/2 



[EK — v)n\ 



G{r) 



fi(K — u)T (n + 2z/) 
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2u+l 



fl f2 



X ( I 1 . (2.4) 

Here, (^) are the Laguerre polynomials given by flA.2p and we use the following notations: 

K = ±{j + 1/2) , u = v/k2 -/i2, ^ = aZ = Ze^/hc, (2.5) 
a = Vl — e = E/mc^, P = mc/h=l/\, 

and 

^ = 2apr = 2Vl^ ^ r. (2.6) 
The elements of 2 x 2-transition matrix in ( 12. 4p are given by 

. a/^ . a (k - z/) 

Ji = , 12 = fi-i^, 91 = , 92 = At- (2.7) 

en — u SK — u 

This particular form of the relativistic radial functions is due to Nikiforov and Uvarov [ID]; it is 
very convenient for taking the nonrelativistic limit c — ^ oo (see also [60j). 

The relativistic discrete energy levels e = En = En/ Eq with the rest mass energy Eq = mc^ are 
given by the Sommerfeld-Dirac fine structure formula 

En = , = . (2.8) 



I • 



Here, n = rir = 0,1, 2, ... is the radial quantum number and k = ± (j ' + 1/2) = ±1, ±2, ±3, ... . 
The following identities 

en = aiy + n) , efi + au = a {n + 2z/) , en — au = an, (2.9) 

2 2 2 2/io\ 2 2 2 

en — V = a n[n + lu) = jj, —an 
are useful in calculation of the matrix elements below. 

The familiar recurrence relations for the Laguerre polynomials allow to present the radial functions 
([221) in a traditional form [2], [7], [I5], [M], [60] as follows 



Fir) 



a 



2, 



IS'^W— e'^-ie-^/^ (2.10) 



G{r) j \j n{K-iy) {eK -i^)T{n + 2p) 

/3i P2)\ Lf iO 



X 



where 

cti = Vl + e (^(k - z/) Vl + + /iVl - , a2 = -\/l + e ({k - u) + e - /iVl - , (2.11) 

p^ = ,/l^(^l^i^-u),/l^ + ^,/TZ^Y /52 = v^T^ ((/t - z^) - yuv^r^) (2.12) 

and a convenient identity holds 

(^{k - u) VTTe ± /iVl - ^ = 2 (k - z/) (k - z/e ± a^u) . (2.13) 
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We give the explicit form of the radial wave functions (12.41) for the lsi/2 state, when n = rir = 0, 
/ = 0, J = 1/2, and k = -1 : 

Here, z/i = a/1 — = ei, = 2^/1 — ef/3r = 2Z (r/ao) , and oq = h^/me^ is the Bohr radius. 
One can see also [2], [7], [8], [16], [19], [25], [28], [31], [37], and [50] and references therein for more 
information on the relativistic Coulomb problem. 



3. Evaluation of the Matrix Elements 
We evaluate the following integrals of the radial functions: 



POO 

Ap = rP+^F^r) + G^r)) dr, (3.1) 
Jo 

POO 

Bp = rP-^^F^r) - G^r)) dr, (3.2) 
Jo 

POO 

Cp = rP+^F{r)G{r) dr (3.3) 



Jo 

in terms of generalized hypergeometric series. (Their relations with the expectation values of the 
operators (Or^), where O = {1, /3, mn/3} , are discussed in the appendix B.) The final results with 
the notations from the previous section can be presented in two different closed forms. Use of the 
traditional radial functions (12.1 Op results in: 

T (2iy + 1) f 1 ~ n, -n, » + 1 \ 

r(2i/ + p + l) " ^ ^ 2iy + l, 2 J ^ ^ 

1 - n, -p, p + 1 \ , ^ f -n, -p, p + 1 

+ (/. + «.) sF,\ ^^^^^ ^ j+ (^-a.) 3F2(^ 2^^^^ ^ 

2, i2apr J^^'^^^^ Bp = 2pan | (3.5) 

T{2u + p + l) " ^ ^ 2^ 2Z/+1, 2 ' ^ ^ 

1-n, -p, p + l\ ( -n, -p, p+1 

+e (u + aK) ■il'o \ + e [u — an) 3_^2 

^ ' 2Z/+1, 1 ) ' \ 2u + l, I 

4/" (2«/3)' Cp (3.6) 



r(2z/ + l) 
V{2u + p+l) 



l-n, -p, p+l\ ^ I -n, -p, p + 1 

ai^i + an) 3F2\ \ - a - an) 3F2 \ 

2z/ + 1, 1 / \ 2z/ + 1, 1 

Nikiforov and Uvarov's form (12.41) gives the following result: 

4/XI/2 {2a(3Y Ap (3.7) 
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, ,r(2z/ + p + 3) / 1 -n, p + 2, -p- 1 

= an [en + v) ; r — ■iho 

^ ^ ' r(2z/ + 2) ' M 2Z/ + 2, 1 



-2 (p + 2) a>ra 



r(2z/ + p + 2) / p + 2, 



V{2u+l) ^ M 2i/ + l, 



,r(2z/ + p+l) /-n, p + 2, 
r(2.) ^^^1^ 2., 1 

4/xi/(2a/3)P Bp (3.8) 

, ,r(2z/ + p + 3) / p + 2, 

= a (£K + z/) , ' 3F2 

^ ^ ^ r(2z/ + 2) ' ' \^ 2Z/ + 2, 1 

,r(2z/ + p + l) f-n,p + 2,-p-l 

— a [SK — Z/ ; r 

^ ^ r(2z/) ' ' 1^ 2z/, 1 

8fiu^{2af3y Cp (3.9) 
r(2z/ + p + 3) 

= + r(2. + 2) ^"^^ 

-2(, + 2)a^.n^i^^±^3F. 

^■^ ' r(2z/ + i) 
^ r(2z/ + p+i) 

+a/i (e/t - z/) ^^^^-^ 3F2 

Here, the terminating generalized hypergeometric series 3F2 (1) are related to the Hahn and Cheby- 
shev polynomials of a discrete variable jSSj, I60j. (See Eqs. (13.111) and flA.ip below, we usually omit 
the argument of the hypergeometric series 3F2 if it is equal to 1.) Two more forms occur if one takes 
one of the radial wave functions from (12. 4p and another one from (12.101) . We leave the details to 
the reader. 

The averages of for the relativistic hydrogen atom were evaluated by Davis [IS] in a form 
which is slightly different from our equations (13.41) and (13. 7p : see also [3] and Ref. [60] for a simple 
proof of the second formula including evaluation of the corresponding integral of the product of two 
Laguerre polynomials: 

^-x^a+s j^a ^3^^g^ 

, T{a + s + l)T{P + m + l)T{s + l) 




3-^2 



ml {n — m)\ F (/? + 1) F (s — + m + 1) 
— m, s + 1, [3 — a — s 
/? + l, n — m + 1 



n> m. 



(The limit c ^ 00 of the integral Ap is discussed in [60j.) Equations (l3.5l) -( l3T6l) and (l3.8l) -( l3T9l) . 
which we have not been able to find in the available literature, can be derived in a similar fashion. 
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It does not appear to have been noticed that the corresponding 3F2 functions can be expressed in 
terms of Hahn polynomials: 

h^:^ (x, N) 



-1 



-X 



n\V{N - n) 



3-^2 



(3.11) 



The ease of handling of these matrix elements for the discrete levels is greatly increased if use is 
made of the known properties of these polynomials [20], [32], and 



For example, the difference-differentiation formulas (4.34)-(4.35) of Ref. [60] (see also (lA.Sp below) 
take the following convenient form 




p{p+l) 



3-^2 




n, 1 — p, p + 2 
2z/ + 2, 2 



(3.12) 



in terms of the generalized hypergeometric functions. (Another proof of these identities is given in 
the appendix A.) As a result, the linear relation holds [51], [1] 



2k {Ap - eBp) -(p + l) {Bp - eAp) = AfxC. 
and we can rewrite fl3.4p - fl3.5p in the following matrix form 



2 (p + 1) afi {2apy 



r (2z/ + 1 




(3.13) 



(3.14) 



where 



■Ji = {fi + an) (a {2eK + p + 1) — 2efi) 
= + an) (a {2k + e {p + 1)) - 2/i) 



72 = (/i — ahi) {a {2eK + p + 1) + 2en) , 
^2 = (yu - an) {a {2k + e {p + 1)) + 2/i) 



(3.15) 
(3.16) 



This representation of integrals Ap and Bp involves the Chebyshev polynomials of a discrete variable 
{x, — 2z/) at a; = n, n — 1 only; see also equation (13.60 for Cp. The corresponding dual Hahn 
polynomials [39] may be considered as difference analogs of the Laguerre polynomials in equation 
(I2.10p for the relativistic radial functions. 



4. Inversion Formulas 

Due to the symmetry of the hypergeometric functions in (13. 40 - 03. 6p under the transformation 
p — > —p — 1, one gets 

^ r(2z/+p+l) {l-6^)p ' ^ ^ 
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~ ^^""^^ T{2u + p+l) {T^^ ' ^ 



= (2a/3)^-+^ ^^^^ Cp. (4.3) 



^2p+i r (2z/ - p) 
r(2z/ + j9 + 1) 

(These relations allow us to evaluate all the convergent integrals with p < —2.) Indeed, 

^-P-i - eB-P-i = (2«/?)'"^' r (2^1^ + 1) " "^^^ ' ^^-^^ 

which gives the first two equations, if Bp ^ eAp and p ^ 0,-1. The last one follows from (13. 6p . 
Special cases p = 0,-1 of (I4.4l) - (l4.5p are simply identity (I6.16P and Fock's virial theorem (I6.13p . 
respectively. In view of our formulas (I3.4I) -( |331) . equation Bp = eAp occurs only when p = or 
n = 0. 

The symmetry of the hypergeometric functions in (I3.7p - (l3.9p under another reflection p —p — 3 
gives 

A_p_, = i2apf^-'' ^('^^-P-^l (4.6) 
^ ^ ^ r(2z/ + p + 3) ^ ^ 

f Af^^ (2p + 3) + (p + 2) {Au^ + {p+l){p + 2)) 

X I ^T) 

V P + 2 

- 2k {2p + 3) Bp- SKfx-^^ Cp\ , 

r (2z/ + p + 3) 
X [-2k {2p + 3) Ap+ (4z/2 + (p + 1) (p + 2)) Bp 



+ 4/i(2p + 3) 



2p+3 r(2z/-p-2) 



^ '^^ r(2z/ + p + 3) ^ 



X ( 2ft:/i ^ ^ ^ - /i (2p + 3) 5^ 



2p + 3 
p + 2 

4/i2 (2p + 3) + (p + 1) (4z/2 - (p + 2)') 



pT2 

as a result of elementary matrix multiplications. These relations can be used for all the convergent 
integrals with p < —3. Further details are left to the reader. 



The corresponding single two-term nonrelativistic relation was found by Pasternack [H], [12] (see 
also pi] and references therein). We have been unable to find the relativistic matrix identities 
( H?T|) - (l4T3l) and (1TO - (1T81) in the available hterature (see Eq. (18) of Ref. [3] as the closest analog). 
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5. Recurrence Relations 

A set of useful recurrence relations between the relativistic matrix elements was derived by 
Shabaev [51j (see also [18], [62j, ^56j, and [1]) on the basis of hypervirial theorem: 

2kAp -{p+l)Bp = AfiCp + 4/3eCp+i, (5.1) 

2KBp-{p+l)Ap = ApCp+i, (5.2) 

fiBp~{p+l)Cp = (3{Ap+i-eBp+i). (5.3) 
Linear relation (13.131) and convenient recurrence formulas 

, 4u'e + 2K{p + 2) + e{p+l){2^e + p + 2) ^ 

^-^^ = -(^ + ') 4(l-e^)(p + 2)/3/. ^'-'^ 



V {p + 2) + {p+ 1) {2k6 + P + 1) {2k6 + P + 2) ^ 

B 



4{l~e^){p + 2)(3fi 



pi 



R , , ,, 4z/^ + 2/.5(2p + 3)+e^(p + l)(p + 2) ^ 

^-^^ = + 4(l-.2)(p + 2)/3/. ^'-'^ 

4^^£ (p + 2) + (p + 1) {2Ke + p + l){2K + e{p + 2)) 
+ 4(i_e2)(p + 2)/5/i ^' 

Cp+i = ^ (2/€ + £ (p + 2)) Ap+i - i- (2/€£ + p + 2) (5.6) 

are obtained from these equations (see [51], [56], and [1] for more details). Their connections with 
the theory of generalized hypergeometric functions will be discussed elsewhere. 

6. Special Expectation Values and Their Applications 

The Sommerfeld-Dirac formula (12. 8p is derived for a point charge atomic nucleus with infinite 
mass and no internal structure (electron moving in static Coulomb field). In reality, the electron's 
mass is not negligibly small compared with the nuclear mass and one has to consider the effect 
of nuclear motion on the energy levels. Actual nuclei have a finite size and possess some internal 
structure, such as an internal angular momentum or spin, a magnetic dipole moment, and a small 
electric quadrupole moment associated with the spin, which also affect the energy levels. Radiative 
corrections are introduced by the quantization of the electromagnetic radiation field. (See [8], [TT] . 
[52] , [53] . [51], [55], [63], and |j58j and references therein for more details.) Calculations of the real 
energy levels of the high-Z one-electron systems with the help of the perturbation theory require 
special relativistic matrix elements. 

From the explicit expressions (13.4p - (l3.9p one can derive the following special matrix elements: 

, 2\ 5n (n + 2i^) + 4z/2 + 1 - £K (2£ft + 3) ,^ ^ , 

^2 = {r) = ——2 (6.1) 

2 (a/5) 

2n'^e^ + 3Ke3 + (S/i^ -p"^ -i)e^ - 3/te - z/^ + 1 
2/5^(1 -£2)2 ' 

2/?/z(l-£2) ' ^^-^^ 
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= (1) = 1, (6.3) 

1\ /5 



^-1 = (j = |;(i-^')(^- + ^^'i^) (6.4) 



/iZ/ 

m^c^ - E"^ ( E I m'^c^-E^ 
+ 



/1\ 2a'(3'K {2eK - 1) 

= \;^/ = 2 ("^^ .(.^-1) (4.^-1) • 

(Note that A^s exists only if |k| > 2 ^Slj-) The average distance between the electron and the 
nucleus f = (r) is given by Ai. The mean square deviation of the nucleus-electron separation is 

{r-ry = A2-(Ai)'. The energy eigenvalue (E) , mean radius (r) and mean square radius (r^) are 
frequently used when making comparisons of wave functions computed by different approximation 
methods. The integrals Ai and A2 have been evaluated in [23], [13], [H], and [60j| (see also Ref. [3] 
for closed-form expressions for {^p}p=_g)- Matrix element A_^ appears in calculation of the electric 
quadrupole hyperfine splitting [13], [52], and [56]. Integrals Ap are also part of the expression for 
the effective electrostatic potential for the relativistic hydrogenlike atom [60] . 

B2 = {pr'^) = — ^ (5n (n + 2z/) + 2z/2 + 1 - 3e/t) (6.7) 
2 (a/3) 

3Ke^ + (5/^2 + 3zy2 _ 1) ^2 _ 3^^ - 3z/2 + 1 

' ^' ' 2/?/i(l-£2) ' ^^-^^ 

i?o = (/?)=£ = (6.9) 

/,\ 

\ r / ze'^mc^ 

^.../4\. ^w-f . (6.11, 

The integral Bq appears in the virial theorem for the Dirac equation in a Coulomb field, 

E = mc^ {(3) , (6.13) 

established by Fock [22] and then developed by many authors (see [H], [12], [15], [32], [19], [T7] . 
[36J, [S], [I8j, [48j, [21], [21], [H], and ^ and references therein). Relation fICT]) can also be 
obtained with the help of the Hellmann-Feynman theorem. 
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(see [18], [36], [5], and [6] and references therein), if applied to the mass parameter [T], [56]. This 
theorem imphes two more relations 

|^..«.C.. ,6.15, 

The following identities hold 

- = ^ = ^ (/^5-2 + (6.16) 

by 05.31) . The integral B_i is evaluated in [TT] and yl_27 -B_25 C~2, and A_3 are given in [JT] 
(see also [SH])- 

The relativistic recoil corrections to the energy levels, when nuclear motion is taken into con- 
sideration, require matrix elements A_2, -B-i and C_2 (see [11], [53], [51], and [1] and references 
therein) . 

na"^ {3n (n + 2z/) + 21/^ + 1) - 3iJ,^e 
4/i (a/5)^ 

K (1 - e^) (1 - z/^) + Sefi^ jeK - 1) 
4/i/?2 (1 - £2) 

= = (2«:E - mc') , (6.18) 

C-i = -^^a^,<J=^C„ (6.20) 
2/iz/ 



*-^2 = . , ^,2 (6-17) 



_ g^/j^ (25/. - 1) _ AjapfC, 

r , ^,3 ^(l-^^)(l-^^)+3e^M^^-l) 4(a/?)^C2 ..... 

/iz/ [y^ — 1) (4z/"^ — 1) V [y^ — 1) (4z/^ — 1) 

The integrals Co, Ci, and -B_i are computed in [2^. In view of (16.101) and f l6.19p . respectively (16.51) 
and (I6.2ip . the following simple relations hold 

Co = ^5-1, A_2 = ^C_2 = ^^^C^C, (6.23) 
" 2/5 ' /i ^ /ii^(4z/2-l) ^ ' 

The last but one was originally found in [12] . 

The integral C\ occurs in calculations of the bound-electron g factor (the anomalous Zeeman 
effect in the presence of an external homogeneous static magnetic field) [33], [17], [6l], [56], and 
[57] . The matrix element C_i has also been found by the Hellmann-Feynman theorem (16.151) . The 
integral C_2 appears in calculation of the magnetic dipole hyperfine splitting [lOj, ^12j, [35j, [21j . 
[13], [S], and [S2]. 
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The author hopes that the rest of matrix elements will also be useful in the current theory of 
hydrogenlike heavy ions and other exotic relativistic Coulomb systems. Professor Shabaev kindly 
pointed out that the formulas derived in this paper can be used in calculations with hydrogenlike 
wave functions where a high precision is required. 



In Table 1, we list the expectation values for the lsi/2 state, when n = Ur - 
and K = — 1. The corresponding radial wave functions are given by Eq. fl2.14p . 

Table 1. Expectation values for the lsi/2 state. 



0, / = 0, i = 1/2, 



P 



"1 2 

-(-'\ z/i K + 1) (2z/i + 1) 



2 \ZJ 



vi + 1) (2z/i + 1) 



(2z/i + 1) 



_2Z 



z/i(2z/i + i: 



2an 



1^1 



Z 
ao 



X 

2ul 



-2 



yi (2z/i - 1: 



(2z/i - 1] 



' Z' 



A 



zy 2 



zy 2 



ao 



A 



z/i (z/i - 1) (2z/i - 1) 



In the table, = z/i = ■>/! — yU^ = y 1 — « = e^/^c is the Sommerfeld fine structure 

constant, ao = /me^ is the Bohr radius, and A = h/mc is the Compton wavelength. The relations 

\Z 



B„ 



El Ap, 



Cn 



2an 



A 



pi 



^ \2Z) r(2z/i + l) 



(6.24) 



(for all the suitable integers p > —2ui — 1 > —3) follow directly from (13.41) . (13. 5p and (13.71) . (13. 9p . 
(The formal expressions for A_3, -B_3, and C_3, when the integrals diverge, are included into the 
table for "completeness"; see Ref. [I] for more details.) The reflection relation (14. 3 p holds for all 
the convergent integrals Ap, Bp, and Cp. 
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Appendix A. Generalized Hypergeometric Series 



, bq-, Z) 



The generalized hypergeometric series is defined as follows [201 

ai' ^2, ... , ap ^ ^ ^ ^ ^ (^ai)^ {a2)„ ... (a 



pFq (ai, 02, 



pJ- q 



b ' 



n=P 



(A.l) 
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where (a)^ = a {a + 1) ... {a + n — 1) = F (a + ra) /F (a) . In this paper we always have p = 3, q = 2, 
z = 1, and tti is a negative integer when the series terminates. The Laguerre polynomials are given 

by m, m, m- 



r(a + n + r 



-n 



The required identity fl3.12p can be derived from the theory of classical polynomials in the follow- 
ing fashion. Let us start from the difference equation for the Hahn polynomials i/m = hm' (x, N) 
[39]: 

{cr{x)V + T{x))Ay^ + X^U^ = 0, (A.3) 
where A/ (x) = V/ (x + 1) = / (x + 1) — / (x) and 

a (x) = X (a + - x) , (A.4) 
r(x) = {P + 1){N -l)~{a + p + 2)x, 
Am = m (a + /3 + m + 1) , 

and use the familiar difference-differentiation formula: 

A/i(^' l'){x,N) = {a + P + m + 1) "-"'^ (x, AT - 1) . (A.5) 

As a result, 

(a (x) V + r (x)) • (x, AT - 1) + m/i^^' (x, A^) = 0. (A.6) 
Letting a = P and /3 — > — 1, one gets 

X (A^ - X - 1) V/iJ^L? (x, A^ - 1) = -m lim h^^' ^'^ (x, A^) (A.7) 



(-l)'"m(m-l)-^- ^x 3F2 



r (A^ - 1) / 1 - m, m, 1 - X 

F(Ar-m)^ '^M 2, 2 - AT 



by (13.111) . The last identity takes the form (13.121) . if the Chebyshev poljTiomials of a discrete variable 
h^m-i (^) A^ — 1) are replaced by the corresponding generalized hypergeometric functions. (Use of 
(lA.Sp in (1A.7P gives the special 3F2 transformation.) 



Appendix B. Dirac Matrices and Inner Product 
We use the standard representations of the Dirac and Pauli matrices: 



with 

0/' \ 1 



0H°° , 1- (B.3) 
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The inner product of two Dirac (bispinor) wave functions 



13 



Ui 
Vl 



^3 



U2 
V2 



V2 
03 
V04/ 



(B.4) 



is defined as a scalar quantity 



R3 



tp'^cj) dv = J ^u|u2 + vjv2 j dv 

(^101 + ^202 + ^303 + V^404) 



(B.5) 



and the corresponding expectation values of a matrix operator A are given by 

{A) = {i;, A^j) . 

From this definition one gets 

where the integrals Ap, Bp, and Cp are given by fl3.1l) - (l3.3p . respectively. 



(B.6) 
(B.7) 



Indeed, the first relation is derived, for example, in Ref. [60] and the second one can be obtained 
by integrating the identity 



r' X*) ( J _°i ) ( ^ ) = (¥■*. X*) ( 4 ) (B.8) 

(we leave details to the reader) in a similar fashion. 



In the last case, we start from the matrix identity 
(an) (3il) = 
and use the Ansatz [60j 



crn 
crn 



-X 



(o-n) cp 



(B.9) 



cp = cp (r) = y (n) F (r) 



-?((^Tn)3^(n)) G(r) 



(B.IO) 



where n = r/r and 3^ = (n) are the spinor spherical harmonics given by (12.31) . As a result. 



irV((an) piP) 



(B.ll) 



{<fi\ x') 



-(o-n) X 
(crn) (p 



ir'P {Fy\ iGy^ (o-n)) ' '^'^ 



(crn) 



= -rP {y^y) FG - rP {y^ (o-n)' 3^) FG = -2r^ (y^y) FG 

with the help of the familiar identity (crn)^ = n^ = 1. Integration over in the spherical coordi- 
nates completes the proof. 
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